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Anomalous periodicity of the current-phase relationship of grain-boundary Josephson junctions 

in high-7 c superconductors 

E. IFichev, V. Zakosarenko, R. P. J. Usselsteijn, H. E. Hoenig, V. Schultze, and H.-G. Meyer 
Department of Cry oelectronics, Institute for Physical High Technology-, P.O. Box J 00239. D-07702 Jena, Germany 

M. Grajcar and R. Hlubina 
Department of Solid State Physics, Comenius University, Mlynska Dolina F2 t 842 15 Bratislava, Slovakia 

(Received 28 January 1999) 

The current-phase relation (CPR) for asymmetric 45° Josephson junctions between two rf-wave supercon- 
ductors has been predicted to exhibit an anomalous periodicity. We have used the single-junction interferom- 
eter to investigate the CPR for these kinds of junctions in YBa 2 Cu 3 0 7 _ x thin films. A remarkable amplitude of 
the 7r-periodical component of the CPR has been experimentally found, providing an additional source of 
evidence for the d-wave symmetry of the pairing state of the cuprates. [SO 163- 1829(99)05629-5] 



A number of experimental results confirm d x t- y i -wave 
symmetry of the pairing state of high-temperature 
superconductors. 1 An unconventional pairing state requires 
the existence of zeros of the order parameter in certain di- 
rections in momentum space. Thermodynamic and spectro- 
scopic measurements do indeed suggest their existence, but 
by themselves they do not exclude conventional s-wave pair- 
ing with nodes. 1 Direct evidence for the d- wave pairing state 
is provided by phase-sensitive experiments, which are based 
on the Josephson effect. 2 Quite generally, the current-phase 
relationship (CPR) of a Josephson junction, I(<p) is an odd 
periodic function of <p with a period 27r. 3 . Therefore I(<p) 
can be expanded in a Fourier series 



I(<p) = I\ sin<p + / 2 sin2<p+- 



(1) 



In the tunnel limit we can restrict ourselves to the first two 
terms in Eq. (1). Since the order parameter is bound to the 
crystal lattice, f(<p) of a weak link depends on the orientation 
of the t/-wave electrodes with respect to their boundary. The 
existing phase-sensitive experiments exploit possible sign 
changes of I\ between different geometries. 2 In this work we 
present a phase-sensitive experimental test of the pairing 
state symmetry of cuprates. Namely, in certain geometries, 
the I x term should vanish by symmetry. In such cases, the 
CPR should exhibit an anomalous periodicity. 

Let us analyze the angular dependence of I x2 in a junction 
between two macroscopically tetragonal rf-wave supercon- 
ductors. As emphasized in Ref. 4, also heavily twinned 
orthorhombic materials such as YBa^^Oy-j belong to this 
class, if the twin boundaries have odd symmetry. We con- 
sider an ideally flat interface between two superconducting 
electrodes. Let 6 X (0 2 ) denote the angle between the normal 
to the grain boundary arid the a axis in electrode 1 (2), see 
Fig. 1 . If we only keep the lowest-order angular harmonics, 
the symmetry of the problem dictates that 4 



/i=/ r cos20| cos20 2 + Asin20i sin20 2 . 



(2) 



The coefficients I c J s are functions of the barrier strength, 
temperature T, etc. The I 2 term results from higher-order 
tunneling processes and we neglect its weak angular depen- 



dence. It is seen from Eq. (2) that the criterion for the obser- 
vation of an anomalous period of the CPR, /j =0, is realized 
for an asymmetric 45° junction, i.e., a junction with & x 
=45° and 0 2 = O. 

The I 2 term is also present in weak links based on con- 
ventional 5-wave superconductors but for all known types of 
weak links |/ 2 //i|<l. For instance, for a tunnel junction 
|/ 2 //)|<n. For a superconductor-normal-metal - 
superconductor (SNS) junction, /asin<p/2 at T=0* and the 
Fourier expansion of Eq. (1) leads to 7 2 //i = -2/5. There- 
fore a possible experimental observation of |/ 2 II x \> 1 in an 
asymmetric 45° junction provides direct evidence of */-wave 
symmetry of the pairing state in the cuprates. 

We have investigated the CPR of YBa 2 Cu 3 0 7 _ x thin-film 
bicrystals with asymmetric 45° [001]-tilt grain boundaries as 
sketched in Fig. 1, using a single-junction interferometer 
configuration in which the Josephson junction is inserted into 
a superconducting loop with a small inductance L. In a sta- 
tionary state without fluctuations, the phase difference <p 
across the junction is controlled by applying an external 
magnetic flux <t> e penetrating the loop: <p= <p € - Pf(<p). Here 
<p e =27T<$> e /<& 0 ; 4> 0 =2.07X 10" 15 Tm 2 is the flux quantum; 
/(<p) = /(<p)// 0 is the CPR normalized to the maximal Jo- 
sephson current / 0 , and £=27rZ,/ 0 /<I> 0 is the normalized 




FIG. 1. Washer-shaped interferometer with one short Josephson 
junction (not in scale). Dimensions are given in the text. 
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critical current. In order to obtain the CPR for the complete 
phase range — 7r^<p^7r the condition /?<! has to be ful- 
filled, because for /?> 1 the curve <p((p e ) becomes multival- 
ued. Following Ref. 3, we express the effective inductance of 
the interferometer using the derivative /' with respect to <p 
as L inl =L[ \ + !/#/"'(<?)]. The inductance can be probed by 
coupling the interferometer to a tank circuit with inductance 
L J-, quality factor Q, and resonance frequency o> 0 through 
the mutual inductance M. 8 External flux in the interferometer 
. is produced by a current I dc + l r f m me ten k co *l am * can De 
expressed as <p e =2TT(l dc + I rf )M/Q> 0 =<p dc +<p rf , where 
M 2 = k 2 LL T with k a coupling coefficient. Talcing into ac- 
count the quasiparticle current in the presence of a voltage V 
across the junction the phase difference is given by the rela- 
tion <p=<Pdc+<Prf-Pf(<P)-2'nT(<P)W®o, where ^^P) 
-L(Rj{ip) with Rj(<p) the resistance of the junction. In the 
small-signal limit <p r /<l and in the adiabatic case wr^l, 
keeping only the first-order terms, the effective inductance 
L € ff of the tank circuit-interferometer system is 



k 2 /3f(<p) 



L eff =L T \\-k* — \=L T \ 1- 



Thus the phase angle a between the driving current and the 
tank voltage U at the resonance frequency of the tank circuit 
o) 0 is 



tan a ((f) = 



k 2 Q0f'(<p) 
1+0T(<P)' 



(3) 



Using the relation [1 +#/"'( <p)] d<p-d<p dc which is valid for 
(p r /<\ and a)T<€\, one can find the CPR from Eq. (3) by 
numerical integration. 

The advantage of the CPR measurement of an asymmetric 
45° junction with respect to the by-now standard phase- 
sensitive tests of pairing symmetry based on the angular de- 
pendence of /] is twofold. First, it avoids the complications 
of the analysis of experiments caused by the presence of the 
term /, . 4 Second, flux trapped in the interferometer washer 
(see Fig. 1) does not invalidate the conclusions about the 
ratio \I 2 H\ | and hence about the pairing symmetry, which is 
not the case in standard phase-sensitive tests of the </-wave 
symmetry of the pairing state. 9 

The films of 100-nm thickness were fabricated using stan- 
dard pulsed laser deposition on (001) oriented SrTi0 3 bi- 
crystalline substrates with asymmetric [001] tilt misorienta- 
tion angles of 45° ±1°. The films were subsequently pat- 
terned by Ar ion-beam etching into 4x4-mm 2 square 
washer single-junction interferometer structures (Fig. 1). The 
widths of the junctions were 1-2 /*m. The square washer 
holes had a side lengths of 50 /im. This geometry of the 
interferometer gives 1^80 pW. The resistance of a similar 
single junction (without interferometer loop) was measured 
directly and R j> 1 ft was found. Therefore the condition 
for the adiabatic limit <or<l is satisfied. For measurements 
of a(<p dc ), several tank circuits with inductances 0.2-0.8 /xH 
and resonance frequencies 16-35 MHz have been used. The 
unloaded quality factor of the tank circuits 10<Q< 150 has 
been measured at various temperatures. The coupling factor 
A* was determined from the period AF dc of a(I dc ) using 
A/A/ rfc = 4>o. Its value varied between 0.03 and 0.09. The 
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FIG. 2. Left panel: Phase angle between the driving current and 
the output voltage measured for sample No. 1 at different tempera- 
tures as a function of the dc current I dc . The curves are shifted 
along the y axis and the data for T- 30 and 40 K are multiplied by 
factor 4 for clarity. From top to bottom, the data correspond to T 
= 4.2, 10, 15, 20, 30, and 40 K. The data measured for 36° bicrys- 
tals (& l t *>36 o ,0 2 t **O) at f=40 K. in the same washer geometry are 
shown for comparison (open circles). Right panel: The same for 
sample No. 3. From top to bottom, the data correspond to T=4.2, 
10, 15, and 20 K. 

amplitude of I r j was set to produce a flux in the interferom- 
eter smaller than 0. 1 4> 0 to ensure the small-signal limit. 

The measurements have been performed in a gas-flow 
cryostat with a five-layer magnetic shielding in the tempera- 
ture range 4.2^ 7"<90 K. The experimental setup was cali- 
brated by measuring interferometers of the same size with 
24° and 36° grain boundaries. We have studied six samples, 
out of which for four samples the 7r-periodic component of 
I((p) was experimentally observed. At low temperatures for 
two samples (Nos. 1 and 2) the value of I 2 is larger than /,. 
For sample Nos. 3 and 4 U is approximately 10-20% of l\ 
and for sample Nos. 5 and 6 1 2 is negligible. As an example 
we plot the phase angle a as a function of the dc current I dc 
for sample Nos. 1 and 3 (Fig. 2). The behavior of sample No. 
1 at low temperatures is defined by the it periodic compo- 
nent of I(<p). The curves for sample No. 3 are 2 7r-periodic, 
nevertheless for the curve at T=4.2 K the local minima 
clearly show the presence of a 7r-periodic component. 

In order to determine the CPR we assume that the period 
of a(I dc ) at T=40 K and A/ rfc = 9.6 /-tA, corresponds to 
A<p rfc = 27r. We take <p dc -0 at a maximum or minimum of 
a. This is necessary in order to satisfy 7(<p = 0) = 0, as re- 
quired by general principles. 3 The experimentally observed 
shift of the first extreme of a(f dc ) from I dc = 0 (Fig. 2) can 
be due to flux trapped in the interferometer washer. Most 
probably, this flux resides in the long junction originated by 
the grain boundary crossing the washer of the interferometer. 
This long junction does not play an active role because the 
Josephson penetration depth is much smaller than the junc- 
tion length, and external fields produced by I dc are smaller 
than the first critical field. Nevertheless, the long junction 
sets the phase difference for I dc = 0 at the small junction. 

In Fig. 3, we show the CPR determined from the data in 
Fig. 2. For all curves we have performed a minimal neces- 
sary shift consistent with 7(<p=0) = 0. Thus we have as- 
sumed that at (p dc = 0 a minimum of a(<p dc ) is realized. For 
an interferometer with a conventional s-wave weak link (and 
also for the 36° junction), at <p dc =0 one gets a maximum of 
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FIG. 3. (a) Josephson current through the junction for sample 
No. 1 as a function of the phase difference (p, determined from the 
data in Fig. 2. The scattering of a((p) values was reduced by fold- 
ing the data back to the interval (0,-n-) and taking the average. Here, 
the symmetry a(<p)=a(-(p) was assumed, (b) The same for 
sample No. 3. The data for the asymmetric 36° bicrystal at T=40 K 
(open circles) are also shown. 

^(<Pdc)- Note that the minimum of a(<p dc ) at <p dc — 0 implies 
a paramagnetic response of the interferometer in the limit of 
small applied fields. 

The amplitude of the ^-periodic component of the CPR 
decreases drastically with increasing temperature, and at T 
= 40 K its contribution is negligible for all samples. The 
temperature dependence of I x and I 2 could be determinated 
with acceptable accuracy for sample No. 1 only. With de- 
creasing T, |/ 2 | grows monotonically down to 7*= 4.2 K, 
while the I x component exhibits only a weak temperature 
dependence (Fig. 4). 
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FIG. 4. Temperature dependence of the Fourier expansion coef- 
ficients l\j determined from the experimental data in Fig. 3(a). 
Solid lines are the Fourier expansion coefficients for the numerical 
data in Fig. 5. 
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FIG. 5. J((p) calculated according to Eq. (64) of Ref. 11 for a 
junction with ^=45.5°, 0 2 = 0, \</=1.5, k=0.5, and T c =60 K. 
I((p) at T=40 K for the 36° bicrystal (open circles) was calculated 
with the same parameters except for 6 X = 36°. 

Our experimental results can be understood as follows. It 
is well known that the microstructural properties of the grain 
boundaries, especialy 45° boundaries, are defined by their 
faceted nature. Faceting is an intrinsic property of the grain 
boundaries, 6 ' 7 and, due to d-wave symmetry of the order pa- 
rameter, the properties of the junctions strongly depend on 
the particular distribution of the facets. Small deviations 
from the ideal geometry of the asymmetric 45° junction lead 
to a finite value of /]. Thus for nearly ideal junctions 
\I 2 /I\\>\ at 7— >0. The region T-~ T c can be analyzed quite 
generally within the Ginzburg-Landau theory. Let the elec- 
trodes be described by the (macroscopic) order parameters 
A u = |A|e' V| . 2 . Then the phase-dependent part of the energy 
of the junction is E = a[A l ^l + U.c.] + b[(A l A* 2 ) 2 ^H.c.] 
+ ••■ where a,b, . . . depend weakly on T. ]0 Thus for T 
close to T c we estimate /,oc|A| 2 *(r c - T) and / 2 a |A| 4 
*(T C -T) 2 , leading to |/ 2 //||<^l. With increasing devia- 
tions from ideal geometry |/ 2 ^il decreases. For large 
enough deviations, negligible values of |/ 2 | are expected. 
These expectations are qualitatively consistent with the ex- 
perimental data (see also Fig. 4). 

So far, our discussion was based solely on symmetry ar- 
guments. Let us attempt a more quantitative analysis of our 
data now. Two different microscopic pictures of asymmetric 
45° Josephson junctions between d- wave superconductors 
have been considered in the literature. The first picture as- 
sumes a microscopically tetragonal material and an ideally 
flat interface. 10 * 12 Within this picture, only sample No. 1 can 
be analyzed. Sample No. 2 had I 0 (T= 1.5 K)= 10~ 2 /iA. At 
this temperature only the 7r-periodic component of /( <p) was 
observed. At higher temperatures / 0 was not measurable. 
I((p) for sample No. 1 calculated according to the model of 
Ref. 11 is shown in Fig. 5. The experimental data can be 
fitted within a relatively broad range of barrier heights. How- 
ever, if we require the I(<p) relation of the 36° junction to be 
fitted by the same (or smaller) barrier height as for the 45° 
junction, we conclude the barrier of the 45° junction to be 
rather low. 14 The dependence of /(<p) on T requires a choice 
of r c ss =60 K in the non-self-consistent theory of Ref. 11. The 
reduction from the bulk T c = 90 K is probably due to a com- 
bined effect of surface degradation and order-parameter sup- 
pression at the sample surface. The temperature dependence 
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of the ratio of the rr and 2 ir periodic components in /( <p) is 
seen to be in qualitative agreement with experimental data in 
Fig. 3(a). This is explicitly demonstrated in Fig. 4 where we 
compare the experimentally obtained I x 2 with the results of 
the Fourier analysis of the curves in Fig. 5. The divergence 
of 1 1 as r— >0 is an artifact of the ideal junction geometry 
assumed in Ref. 1 1 . If a finite roughness of the interface is 
taken into account, this divergence is cut off and the experi- 
mental data in Fig. 4 do indeed resemble theoretical predic- 
tions for a rough interface. 12 However, the non-self- 
consistent theory of Ref. 11 is unable to explain the 
experimentally observed steep CPR close to the minima of 
the junction energy [see Fig. 3(a)]. 

In a different approach a heavily meandering interface 
with 9{= dj(x) is assumed. Now, the critical current density 
j c (x) is a random function with a typical amplitude 
(\j c (x)\)~~j c . If the average critical current along the junc- 
tion (J c )<j c , a remarkable ^-periodic component is present 
in the CPR. The relation |/ 2 ^il depends on the distribution 
of j c (x) and can be much larger than one for (J c )<j c . 15,16 
This model qualitative explains the obtained results for all 
samples, however for a quantitative comparison the actual 
microscopic distribution j c (x) should be known. Note that 



also within the picture of Refs. 15 and 16 the rf-wave sym- 
metry of the pairing state is crucial, otherwise the condition 
(jcl^Jc i s difficult to satisfy. 

Our present understanding of /( <p) in the asymmetric 45° 
junction is only qualitative. We cannot say whether the re- 
markable amplitude of the 7r-periodic component of /(<p) is 
dominated by the microscopically flat regions, 13 or due to the 
spatial inhomogeneity of the junction. This issue requires 
further study. 

In conclusion, we have measured the magnetic-field re- 
sponse of a single-junction interferometer based on asym- 
metric 45° grain-boundary junctions in YBa 2 Cu 3 0 7 _ T thin 
films. A large ^-periodic component of 7(<p) has been ex- 
perimentally found, which is in agreement with theoretical 
predictions for J Jt 2_ v 2-wave superconductors. Hence our re- 
sults provide an additional source of evidence for the (/-wave 
symmetry of the pairing state in the cuprates. 
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For various configurations of Josephson junctions incorporating superconductors with unconventional 
order parameter symmetry, such as most high-T c cuprates, deviations from the standard sinusoidal 
current-phase dependence have been predicted. To this point, these deviations have never been observed 
experimentally. We have measured the current-phase relation of high-7 c Josephson junctions, namely, 
YBa2Cu 3 07-, thin film bicrystals, comprising symmetric 45° [001] tilt grain boundaries. The current- 
phase relations of all junctions investigated were found to be extremely nonharmonic, in agreement with 
a Jj(2_ y 2-wave dominated symmetry of the order parameter. [S003 1-9007(98)06674-5] 



PACS numbers: 74.50. +r 

The current-phase relation (CPR) f(<p) is a character- 
istic property of any weak link connecting two supercon- 
ductors. It describes the dependence of the Cooper-pair 
current I p on the phase difference <p of the order parame- 
ters of both superconducting electrodes. In a general form 
it is expressed as 

I P = Icf(<p), -1 </(?)<!, (1) 

I c being the critical current of the weak link. It was 
shown by Josephson [1] that for ideal tunnel junctions 
between conventional superconductors the CPR is sinu- 
soidal, i.e., f(<p) = sin(^). This sinusoidal dependence 
has been confirmed experimentally numerous times for 
standard tunnel junctions between conventional supercon- 
ductors [2], 

Recently, it has been revealed that the order parameter 
of most high-7V cuprates is unconventional, dominated 
by a d x 2- y 2 symmetry component [3-5]. Because of the 
sign change of the order parameter associated with this 
symmetry, strong deviations from the standard sinusoidal 
dependence have been predicted for the current-phase 
relations of various configurations of Josephson junctions 
employing such unconventional superconductors [6-9]. 
In particular, nonharmonic and double-periodic current- 
phase relations are expected for junctions oriented nomi- 
nally perpendicular to the (1 10) direction of one or of 
both electrodes, as well as for junctions for which the 
(110) direction of one of the electrodes is aligned with 
the (100) direction of the other, such as for 45° [001] tilt 
grain boundaries. 

These predictions are highly unusual. Therefore, an ex- 
perimental clarification of the CPR for high-7 c junctions 
for which deviations from a standard harmonic behavior 
are expected is desirable. Such experiments will further 



enhance the understanding of the influence of the order 
parameter symmetry on the properties of grain boundaries 
and high-7V Josephson junctions. In addition, they will 
provide valuable information for the design and use 
of Josephson junction-based circuits, of which many 
characteristics directly depend on the CPR. However, to 
our knowledge, such experiments have not been carried 
out. All available data refer to Josephson junctions 
for which nominally sinusoidal current-phase relations 
are expected. The CPR was measured for weak links 
prepared by ion irradiation [10], for step-edge junctions 
[11,12], and also for 24° bicrystal grain boundaries 
[12]. In nearly all of these cases sinusoidal current- 
phase relations were found. Deviations from a sinusoidal 
dependence have been observed only for one step-edge 
junction, measured at 77 K [12]. These deviations can be 
explained by the influence of thermal noise [13]. 

For these reasons we have investigated the CPR of 
YBa2Cu307-jt thin film bicrystals with symmetric 45° 
[001]-tilt grain boundaries, as sketched in Fig. 1(a). For 
these junctions, strong deviations from a sinusoidal CPR 
are anticipated. 

Following a standard approach [14], the CPR was mea- 
sured using a single-junction interferometer configuration 
in which the Josephson junction is part of a superconduct- 
ing loop with a small inductance L. The phase difference 
<p across the junction is controlled by applying an external 
magnetic flux $ c penetrating the loop: 

<p = <p c - fif((f) + <p n + 27rm . (2) 

Here, <p e =? 27r<I> e /<&o is the external flux normalized to 
the flux quantum 4> 0 (= 2.07 X 10" 15 Tm 2 ). The con- 
stant /3 = 27rL/ c /4>o is the normalized critical current, 
<p n is a term accounting for the effective noise, and m is 
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FIG. 1. (a) Schematic representation of a symmetric 45° 
[001]-tilt grain boundary junction in a d x z- y i superconductor. 
The boundary in the superconducting thin film is meandering, 
leading to the occurrence of v facets, (b) Schematic of the 
measurement setup. The Josephson junction is denoted by JJ, 
and C indicates the capacitance of the tank circuit. The other 
denotations are explained in the text. 



an integer. As capacitive contributions to the loop current 
are insignificant at the measurement frequencies, the junc- 
tion capacitance has been neglected. The quasiparticle 
current is also negligible, for reasons discussed below. 

The superconducting loop is inductively coupled to a 
tank circuit with inductance L T [see Fig. 1(b)]. This tank 
circuit is driven with a current /rf at a frequency w and 
a dc current / dc . Thus, <p t can be expressed as a sum 
of a dc and an rf component <p e = <pd C + ^rf. In this 
arrangement, the effective impedance Z e ff (o>) of the loop- 
tank circuit combination is a function of <p e . As shown 
by Rifkin and Deaver [14], the CPR can be obtained 
from this dependence, provided that <prf <c 1. To obtain 
the CPR for the complete phase range 0 < <p ^ 2ir, the 
condition /? < 1 has to be fulfilled in addition. 

To enhance the accuracy of the measurement, we have 
adapted this common approach and retrieved the CPR 
from a measurement of the ^>dc dependence of the phase 
angle a between the drive current / r f and the tank voltage 
V at the resonant frequency of the tank circuit <oq. As 
described in Ref [12], at <o 0 , the a(^ dc ) dependence is 
related to the derivative of the CPR f'(<p) = df{<p)/d<p 
in the following way: 



tan or (<pdc) = 



k 2 QPf'(<p{<P*c)) 
1 + Pf'(<p(<p dc )) 



(3) 



Here k is the coupling factor between the tank induc- 
tance and the interferometer, k 2 = M 2 /{LLt), where M 
is the mutual inductance [Fig. 1(b)]. Using Eq. (3), from 
the measured a{<pdc) dependence f(<p(<p<ic)) is obtained. 
The CPR is restored by integrating f'(<p(<Pd C )) numeri- 
cally, using the d<p{<Pd C )/d<Pdc dependence obtained from 
differentiating Eq. (2) with respect to ^>dc 

The samples investigated consisted of three bicrys- 
talline YBa 2 Cu 3 0 7 -, films with a T c (R = 0) of 88 K. 
The films, with thickness t = 100 nm, were deposited by 
standard pulsed laser deposition on (OOl)-oriented SrTi03 
bicrystalline substrates [15] with symmetric [001]-tilt mis- 
orientation angles of 45° ± 2° and were subsequently 
patterned by Ar ion-beam etching into 8X8 mm 2 or 
5X5 mm 2 square washer single-junction interferome- 
ter structures. The widths of the junctions were b 
2-3 /zm. The washer holes had a side length of 50 ^tm, 
leading to L « 80 pH. 

To minimize the influence of external noise, the 
samples were measured in superconducting and double 
magnetic shielding at a temperature of 4.2 K. The 
condition /? < 1 for the investigated interferometers 
was confirmed experimentally from the character of its 
response versus <p^c [12]. 

YBa2Cu307-jc grain boundaries with a symmetric [001] 
tilt angle of 45° typically have a normal-state interface- 
resistivity p„ > 1 X 10" 8 ft cm 2 , which we also mea- 
sured for boundaries fabricated under identical conditions 
as the junctions used in the present experiments. In the 
configuration used, this p n corresponds to normal-state 
resistances R n > 1 ft. Accordingly, the relaxation time 
of the interferometer t = L/R is short (t <sc l/a>o). and 
hence the quasiparticle current is negligible [12]. 

For the measurements of a(<p^ c ) 9 two tank circuits with 
quality factor Q = 120 and inductance Lj — 0.4 /zH, 
co 0 = 30 MHz, and L T = 0.73 fiH, (o 0 = 23 MHz, re- 
spectively, were employed. The phase angle was recorded 
as a function of / dc after amplification of the tank voltage 
by a high-impedance amplifier. The coupling coefficient 
k was determined from the period of the a(/d C ) depen- 
dence. Values of k = 0.072 and k = 0.054 for the re- 
spective tank circuits are obtained. To ensure the validity 
of the small-signal limit, the measurements were carried 
out with ^ r f < 0.15. 

A typical a(/d C ) dependence is shown in Fig. 2. The 
corresponding CPR, depicted in Fig. 3, is clearly deviat- 
ing from the standard sinusoidal behavior. Samples fabri- 
cated on different substrates and measured with both tank 
coils followed closely the same behavior. It is empha- 
sized that the experimental setup employed and the pro- 
cedure followed are identical to those used to measure the 
current-phase relations of step-edge junctions and thin- 
film bicrystals with a symmetric [001] tilt of 24°. For 
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FIG. 2. Phase angle a between the driving current and the 
output voltage measured at 4.2 K as a function of the dc current 
/ dc , for an YBa2Cu 3 07_, single junction interferometer circuit 
containing a symmetric 45° [001]-tilt grain boundary. 



all of those samples nominally sinusoidal current-phase 
relations were observed, and all deviations of the apparent 
CPR from a sinusoidal one can be attributed to thermal 
noise [12,13]. 

The measured deviations from a sinusoidal de- 
pendence for the current-phase relations of these 
45° bicrystals are startling. It is important to 
note that the effective Josephson penetration depth 
Aj = [*o/(47r/xoOc)A)]^ 2 5 /im is larger than the 
width of the junction b (narrow-junction limit). Here 
A is the London penetration depth. Although several 
mechanisms are known to cause nonsinusoidal current- 
phase relations for narrow junctions fabricated from 
conventional superconductors, all of these mechanisms 
fail to account for the anomalous dependencies presented. 
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FIG. 3. The normalized current through the junction pf(<p) 
as a function of the phase difference <p restored from the 
measured ar(/ dc ) as shown in Fig. 2. For comparison, the 
function 0 sin(^) with ft = 0.8 is plotted as a solid line. 



First, one potential source of such deviations is thermal 
noise. To evaluate its influence we consider a sinusoidal 
CPR and calculate with Eq. (3) the a(<pd c ) dependence, 
assuming a thermally induced Gaussian spread p{(p„). 
With this, the value of tan a(<p& c ) is given by [12] 

tana(^ dc ) = k 2 Q/3 

COS <p(<p c ,<Pn) 



« 1 + /3 COS (p((Pc,<Pn) 



p{<p„)d<p n . 



Using Eq. (4), minor deviations of the CPR from a stan- 
dard sinusoidal behavior are well described quantitatively 
for 24° boundaries measured at 77 K [12]. However, no 
realistic set of /? and <p n exists to account for the large 
deviations of the CPR observed for the 45° boundaries. 

Second, for weak links with a high current density, 
current-induced suppression of the order parameter in 
the electrodes close to the weak link can give rise to 
nonharmonic current-phase relations [16]. It is unrealistic 
that this effect is the cause for the deviations presented 
here, as the intragrain critical current density exceeds 
10 7 A/cm 2 at 4.2 K and is therefore much larger than 
the grain boundary (j c ) < 4000 A/cm 2 at the same 
temperature. 

Third, several additional mechanisms, described in [2], 
lead to deviations from a harmonic CPR. In all these 
cases, the slope of f(<p) at <p — 0 is smaller than at <p = 
7r, which is in contrast to our results. Therefore, these 
mechanisms cannot explain the current-phase relations 
observed either. 

On the other hand, as will be pointed out in the 
following, the measured CPR can be accounted for 
by the unconventional order parameter symmetry of 
YBa2Cu307-jr and by the microstructural properties of 
the grain boundaries, in particular by their faceted nature 
[17,18]. Interestingly, due to the J J :_ > 2-wave character 
of the order parameter, the faceting has a more significant 
influence on the electronic properties of boundaries with 
a misorientation close to 45° than on boundaries with 
considerably smaller misorientation angles [18]. This 
also concerns the CPR, as the 45° boundaries contain a 
higher density of facets that themselves show anomalous 
behavior. Two kinds of such anomalies, to be considered 
here, are described in the literature. 

First, due to the sign difference of the adjacent lobes of 
the rfjc2_ y 2-wave order parameter, many facets are biased 
with an additional it phase shift (ir facets) [17-19]. 
These phase shifts give rise to unconventional junction 
properties, such as a spontaneous generation of magnetic 
flux in the grain boundary junction [19-21]. As described 
in Ref. [21], the local phase difference <p(x) along the 
grain boundary (0 < x < b) can be written as 

<p(x) = + (5) 

where g(x) is a rapidly alternating function accounting 
for the 7r phase shifts and for the spontaneously generated 
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magnetic flux in the junction, and tlt(x) is the remaining 
slowly varying phase difference reflecting the magnetic 
flux in the interferometer loop. For a narrow junction 
(b < AjX t(t is independent of x. With Eq. (5), the time- 
independent sine-Gordon equation, describing the spatial 
dependence of the local phase difference over the junc- 
tion, becomes 

A5^ = ^sin[* + (6) 

The solution f (jc) of this equation, and thus the pattern 
of self-generated flux, depends on tp. The redistribution 
of this flux by a change of if/ is expected to lead to 
remarkable deviations from a harmonic dependence for 
the CPR measured for the entire junction, also if the local 
CPR is nominally sinusoidal. 

Second, it has been proposed [7,8] that the CPR of 
facets formed by the (110) and by the (100) planes 
of the adjacent grains is periodic with tr and thus has 
a double periodicity as compared to the standard case. 
Transmission electron microscopy has revealed that 45° 
[001] tilt grain boundaries in YI^O^O?-* tend to be 
composed for a considerable part of such facets [22]. For 
the whole junction, this leads to an anomalous CPR: 

/ — 7 c isin^r + 7 c2 sin2^, (7) 

by which the observed CPR can be described. 

In summary, the current-phase relations of grain bound- 
ary junctions with a misorientation of 45° were mea- 
sured with a modified Rifkin-Deaver method. The CPRs 
were deduced from measurements of the phase angle be- 
tween the rf drive current and the rf tank voltage. The 
current-phase relations of the Josephson junctions showed 
pronounced deviations from a harmonic behavior, which 
cannot be accounted for by thermal noise or by other stan- 
dard mechanisms, but are attributed to the d x 2- y 2 -wave 
symmetry of the order parameter and the faceting of the 
grain boundaries. 
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The predominant d-wave pairing symmetry in high-temperature superconductors allows for a 
variety of current-phase relations in Josephson junctions, which is to a certain degree fabrication 
controlled In this Letter, we report on direct experimental observations of the effects of a non- 
sinusoidal current-phase dependence in YBCO dc SQUIDs, which agree with the theoretical description 
of the system. 
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It is well established [1] that the wave function of a 
Cooper pair in most cuprate high-temperature supercon- 
ductors (HTS) has a </-wave symmetry. Its qualitative 
distinction from, e.g., the anisotropic .y-wave case is that 
the order parameter changes sign in certain directions, 
which can be interpreted as an intrinsic difference in the 
superconducting phase between the lobes equal to tr. 

The latter leads to a plethora of effects, such as forma- 
tion of Andreev bound states at surfaces and interfaces in 
certain crystallographic orientations [2-4]. The current- 
phase dependence Is(<f>) in Josephson junctions formed by 
dd junctions, as well as by s d junctions comprised of a 
cuprate and a conventional superconductor, depends both 
on the spatial orientation of the d-wave order parameter 
with respect to the interface, and on the quality of the 
latter [5-9]. Time-reversal symmetry can also be sponta- 
neously violated and thus spontaneous currents generated 
[10-12]. Another effect can be doubling of the Josephson 
frequency [6,13,14]. 

In this Letter, we report on experimental observa- 
tions of strong effects of an unconventional current-phase 
relation on the dynamics of two dd junctions integrated 
into a superconducting interference device (SQUID) 
configuration. 

Since I$(<f>) must be a 2?r-periodic odd function, it can 
be expanded in a Fourier series. In most cases, only the 
first two harmonics give a significant contribution to the 
current: 

/5(<£) = 'cSin<£-/?sin2<£. 0) 

In Josephson systems of conventional superconductors, 
the second harmonic will usually be negligible [15] but 
in dd junctions the second harmonic may dominate. If 
/? > /i/2, the equilibrium state is no longer <f> » 0 but 
becomes double degenerate at <f> — ± arccos(/J/2/i I ) — ♦ 
±tt/2. The system can then spontaneously break time- 
reversal symmetry by choosing either state. Spon- 
taneous currents as well as fluxes can be generated in 
this state. The potential will have the shape of a double 
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well, and there are reasons to believe that it will be 
possible to observe quantum coherence in this system. 
The presence of a second harmonic in the current-phase 
relation (CPR) of a dd junction was confirmed by 
IlMchev et al [8]. 

A nonsinusoidal CPR of the junctions will change the 
dynamics of a dc SQUID [16]. Regarding the junctions as 
magnetically small, the supercurrent through the SQUID 
in the presence of an external flux sa 4> 0 ■ (<f> x /2ir) 
can be written as 

hit* 0J-/i|Sin0-/?,sin(20) 

+ I\ 2 sim> + <f> x ) - /? 2 sin2(<£ + <f> x ). (2) 

The critical current through the SQUID is given by the 
usual expression / C (0J a max^ l s (<f>, <f> x ). The time- 
averaged voltage over the SQUID in the resistive regime 
is readily obtained in the resistively shunted junction 
approximatioa By introducing 8[<f>, <f> x ] = <f> 2 ~ 4>\ and 
applying the same method as in [17] with the necessary 
generalizations, we obtain the following for the average 
voltage over the SQUID: 

-,HH,-f)]-, 

(3) 

Here G\ 2 are the normal conductances of the junctions, 
and 

5 + 4>* + 7T 1 ^ * 5 / 2 ) " A + 5/2)] = 0 (4) 

gives the difference, 5, in phase drops across each junc- 
tion. In deriving (3) and (4), we have assumed that the 
inductance L is equally divided between the SQUID 
arms. We have also neglected the spontaneous magnetic 
fluxes in the dd junctions, due to their small amplitude 
[11,18], Though (4) is only explicitly solvable in the limit 
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L — 0, it always yields 5[— <f> f -<f> x ] = -8[<f> f <£ x ]. This 
means that the usual inversion symmetry is retained 

The results of numerical calculations based on (2) and 
(3) are shown in Fig. 1. The cusps in the critical current 
correspond to the points at which the global maximum in 
(1) switches from one local maximum to another [16]. 
Note the quasi-3> 0 /2 periodicity of the current isolines in 
the V - <}> x picture, reflecting the current-phase depen- 
dence (1), and their shift along the <b x axis, which de- 
pends on the sign of the bias current (as it must to 
maintain the central symmetry with respect to the ori- 
gin). The shift does not depend on the magnitude of the 
current since we neglect the self-inductance. For large 
biases, the <E> 0 periodicity is restored Indeed, as the bias 
grows, one set of minima of the washboard potential, 
U = (h/2e)[-I l cos<f> + (/"/2)cos2<£ - disappears 
first unless the first harmonic 7 1 is exactly zero. 

We have fabricated and studied a large number of dc 
SQUIDs. The samples were fabricated from 250 nm thick 
YBCO films deposited on SrTi0 3 bicrystals. The grain- 
boundary junctions (GBJs) are of the asymmetric 
[001]-tilt type with the misorientation angle of 45° 
(0°-45° GBJ). For more information on GBJs see, for 
example, Ref. [19]. 

The pattern was defined using E-beam lithography and 
then transferred to a carbon mask employing a multistep 
process. Finally, the YBCO is etched through the mask 
using ion milling. This scheme allows us to fabricate 
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FIG. 1. The results of simulations of the I c - <f> x and V - <f> x 
dependence for a dc SQUID with I l cl = 1, I l c2 = 0. 1, » 0.2 
and = 0.4 (arbitrary units). The different curves correspond 
to bias currents in the range / = l\ x to / = 5/^,. We assume 
L —0 and G\ = C 2 
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high-quality bicrystal junctions as narrow as 0.2 /im, as 
has been reported elsewhere [20]. In the SQUIDs under 
investigation, the junctions are nominally 2 yum wide; 
hence, the fabrication-induced damage of the junctions 
is small. 

The measurements were done in an EMC-protected 
environment using a magnetically shielded LHe cryostat 
However, the magnetic shielding is imperfect, as is evi- 
dent from the fact that the expected zero-field response of 
our SQUIDs is not exactly at zero. The measuring elec- 
tronics is carefully filtered and battery powered whenever 
possible. In order to measure the dependence of the 
critical current on the applied field, we used a voltage 
discriminator combined with a sample-and-hold circuit 
All measurements reported here were performed at 4.2 K. 

The SQUID loops are (15 X 15) /tm 2 . The numerically 
calculated inductance [21] is approximately 25 pH, yield- 
ing the factor f$ = 27rLI c /® 0 between 0.5-2. 

The SQUIDs were largely nonhysteretic with a resis- 
tance of about 2 fl. The measured critical current varies 
from sample to sample but is in the range of tens of 
microamperes giving a current density of the order of 
J c = 10 3 A/cm 2 . The estimated Josephson penetration 
length kj = &o/y/4irft 0 j c A L ls approximately 2 jxm in 
all junctions, which means that the junctions are magneti- 
cally short This is supported by the quasiperiod of the 
pattern in Fig. 2 being close to the expected value 
<f> 0 /2\ L w [17]. The differential conductance curves do 
not show any trace of a zero bias anomaly (ZBA), as is 
expected for 0°-45° GBJs. ZBAs have been observed by 
other groups in GBJs with other orientations [2], 

The critical current is plotted as a function of applied 
magnetic field for two SQUIDs in Fig. 3. The result is in 
qualitative agreement with theory if we assume that the 
SQUID junctions have different ratios of the first and 
second harmonics of the critical current This assumption 




Applied Magnetic Field (ml) 



FIG. 2 Critical current as a function of magnetic field at 
4.2 K. The dashed box indicates the area plotted in Fig. 3(a). 
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is supported by the fairly small modulation depth [it is 
easy to see from Eq. (2) that I c would go exactly to zero in 
a SQUID with junctions of identical /^/AriL 

We can fit the data to Eq. (2), if we compensate for the 
residual background magnetic field and assume that we 
have a small excess current (of the order of a few jxA) in 
the junctions. The fitting parameters again confirm that 
there is a large asymmetry between the arms of the 
SQUIDs. Note that the model does not consider the flux 
penetration into the junctions, 

The result for fields of the order of mT is presented in 
Fig. 2, which shows the I c modulation of the SQUID 
enveloped by an anomalous Fraunhofer pattern quite 
similar to what has been reported by other groups 
[22,23] for 0°-45° GBJs. Note the inversion symmetry 
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FIG. 3. Critical current as a function of applied magnetic 
field for two different SQUIDs that are nominally identical. 
The solid line represents the fitted expression. The fitting 
parameters are as follows: (a) /J., = 9 fiA, l\ 2 =0.3 ^A, /" = 
3.7 /iA, and 7» 2 = 22.7 /iA; (b) l\ x = 7.8 /iA, l l c2 = 3.0 pA, 
/J. 1 , — 5.3 ftA, and /° 2 = 4.3 /iA. In both cases, the fit has been 
adjusted with respect to the residual background field and the 
excess current of the junctions. 



of the pattern with respect to the origin. That the global 
maximum is not in the center can be explained in several 
ways; it has been shown, for example, that this could be 
due to the presence of so-called 7r loops in the junction 
interface [24]. 

Figure 4 shows the V-B dependence of one of the 
SQUIDs. The pattern is again field inversion symmetric. 
The overall structure is the same as in the model depen- 
dence of Fig. 1, but there is also an additional shift due to 
self-field effects, which depends on the magnitude of the 
bias current and corresponds (at maximum) to a flux 
~-0. l<I>o. In a beautiful experiment, a similar dependence 
was recently observed by Baselmans et al in a Nb-Ag- 
Nb SNS junction where current injectors were used to 
change the occupation of current-carrying states in the 
normal region [25]. A deviation from the model occurs at 
V = 100 £iV where the minima and maxima switch. This 
is probably due to an LC resonance in the SQUID. Taking 
L = 25 pH, this would require C = 0.8 pF, which agrees 
with our measurements on single junctions 

Remarkably, the observed offset of the V-B character- 
istics with respect to the direction of the bias current 
appears to be a much more robust manifestation of the 
presence of a second harmonic of the Josephson current 
than the shape of the I c - B curves itself. We observed the 
shift even in SQUIDs with the smallest junctions down to 
0.5 fiim wide, where the deviations from the usual sinu- 
soidal CPR were not obvious from the I c — B dependence. 

Generally, the nature of the transport through a GBJ 
will depend on its transmissivity D. IHchev et al [8] 
have reported values of D as high as 0.3 in symmet- 
ric (22.5°-22.5°) dd junctions as opposed to the usual 
estimate for a GBJ, D - 10~ 5 -10" 2 . Since usually 
I X c/I l c a A a high-transmissivity GBJ is required in order 
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FIG. 4. Voltage modulation as a function of applied magnetic 
field for the SQUID whose l c - B is shown in Fig. 3(a). The 
pattern is again inversion symmetric. Note the sign change at 
100 /*V, which we believe is due to a LC resonance in the 
SQUID loop. 
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to observe effects of the second harmonic. An estimate of 
the average transmissivity of our junctions would be 
p ab l/R N A ~ 10~ 2 [26] assuming /, the mean-free path, 
to be equal to 10 nm and a resistivity in the a-b plane p ab 
equal to 10" 4 ft cm. This is still too low to explain the 
strong second harmonic we observe. However, it is known 
from, e.g., TEM studies [19], that the grain-boundary is 
far from uniform; the properties can significantly vary 
depending on the local properties of the interface, effects 
such as oxygen diffusion out of the GB, etc., which are 
difficult to control. It is therefore reasonable to assume 
that there are many parallel transport channels through 
the GB [27,28]. Channels with high-transmissivity domi- 
nate the transport and might have D ~ 0. 1 even though 
the average transmissivity is much lower. This is also 
consistent with the fact that most of our SQUIDs seem 
to be highly asymmetrical which is to be expected if the 
distribution of channels is random. The ratios of l\ and /? 
can vary as much as 10 times between two junctions in the 
same SQUID, even though the fluctuations of the total I c 
from sample to sample are much smaller. It is also clear 
from general considerations that a high value of /? ex- 
cludes a high value of since the second harmonic 
usually dominates if the odd harmonics of the super- 
current are canceled by symmetry [29]. 

Recent studies of 0°-45° GBJs have demonstrated that , 
the SQUID dynamics can be altered by the d-wave order 
parameter in YBCO [30]. It is, however, important to 
point out that our results do not directly relate to, e.g., 
tetracrystal 7T-SQUID experiments; the latter crucially 
depend on having one n junction with negative critical 
current, but still only the first harmonic present in I c (<f>). 
Our SQUIDs have a conventional geometry, but uncon- 
ventional current-phase relations. 

One explanation for the pronounced effects of the 
second harmonic could be that relatively large sections 
of the interface are highly transparent and have a low de- 
gree of disorder. This in turn could be related to our fabri- 
cation scheme which seems to preserve the integrity of 
the barrier. This makes feasible their applicability in the 
quantum regime and supports our expectations that quan- 
tum coherence can be observed in this kind of structures. 

In summary, we have observed a very pronounced sec- 
ond harmonic in the current -phase relation of a "conven- 
tional" YBCO dc SQUID with 0°-45° grain-boundary 
junctions. It has strongly influenced the SQUID dynam- 
ics. All details of the SQUID behavior were explained 
within a simple model of a dd junction with relatively 
high transparency. We believe that these effects are im- 
portant for better understanding of HTS Josephson junc- 
tion and SQUIDs. 
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We present various concepts and experimental procedures to produce biepitaxial YBa2Cu307- x 
grain boundary Josephson junctions. The device properties have an interesting phenomenology, 
related in part to the possible influence of "7r-loops n . The performance of our junctions and Su- 
perconducting Quantum Interference Devices indicates significant improvement in the biepitaxial 
technique. Further, we propose methods for fabricating circuits in which "0-" and u 7r-loops" are 
controllably located on the same chip. 



I. INTRODUCTION 

The possibility of realizing electronic circuits in which 
the phase differences of selected Josephson junctions are 
biased by n in equilibrium is quite stimulating. 1 The 
concept of such 7r-phase shifts was originally developed 
in the "extrinsic" case for junctions with ferromagnetic 
barriers 2 and in the "intrinsic" case for junctions exploit- 
ing superconductors with unconventional order parame- 
ter symmetries. 3 As a result of the possible d x 2_ y 2 or- 
der parameter symmetry of high critical temperature su- 
perconductors (HTS), 4 the presence of intrinsic 7r loops 
has also been considered for HTS systems. 5 This has 
been discussed recently in view of novel device con- 
cepts, and in particular for the implementation of a 
solid state qubit 1,6-8 and for Complementary Joseph- 
son junction electronics. 9 In this paper we discuss how 
YBa 2 Cu307-r (YBCO) structures made by the biepi- 
taxial technique 10,11 can be successfully employed to pro- 
duce arbitrary circuit geometries in which both "0" and 
7r-loops are present, and possibly to obtain a doubly de- 
generate state. 1,6 Of course, great caution should be used 
because of stringent requirements on junctions parame- 
ters for practical applications of such devices. 

Josephson junctions based on artificially controlled 
grain boundaries have been widely employed for funda- 
mental studies on the nature of HTS. 4,7,8 The lack of a 



reliable technology based on the traditional trilayer con- 
figuration (i.e. a sandwich type junction with an insu- 
lator between the two superconducting electrodes) also 
enhanced interest in GB Josephson junctions for appli- 
cations. Although the mechanism of high-Tc supercon- 
ductivity and the influence of grain boundaries on the 
transport properties are not completely determined, re- 
producible and good quality devices are routinely fab- 
ricated. YBCO GB junctions are usually classified as 
bicrystals, 12 biepitaxials, 11 and step-edges, 13 depending 
on the fabrication procedure. The bicrystal technique 
typically offers junctions with better performances and 
allows in principle the realization of all different types 
of GBs ranging from [001] and [100] tilt to [100] twist 
boundaries. GB junctions based on the step-edge and 
biepitaxial techniques offer the advantage, with respect 
to the bicrystal technology, of placing the junctions on 
the substrate without imposing any restrictions on the 
geometry. A comparison between the different GB tech- 
niques is far beyond the aim of this paper. Nevertheless 
we intend to show that significant improvements with re- 
spect to the original technique developed by Char et al. 11 
are possible for biepitaxial junctions, and that the result- 
ing devices have potential for applications. As a matter 
of fact, in traditional biepitaxial junctions, the seed layer 
used to modify the YBCO crystal orientation on part of 
the substrate produces an artificial 45° [001] tilt (c-axis 
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tilt) GB. The nature of such a GB seems to be an intrinsic 
limit for some real applications. A convincing explana- 
tion has been given in terms of the d-wave nature of the 
order parameter and more specifically by the presence 
of 7r-loops. 14 As demonstrated by studies on bicrystals, 
based on the same type of 45° [001] tilt GB, the presence 
of 7r -loops reduces the lcR-N values (where Ic and 
are the critical current and the high normal state resis- 
tance respectively), produces a dependence of the critical 
current Ic on the magnetic field H quite different from 
the FVaunhofer-like pattern, and generates unquantized 
flux noise at the grain boundary. 14 

We will show that the implementation of the biepi- 
taxial technique 10 we developed to obtain 45° [100] tilt 
and twist (a-axis tilt and twist) GBs junctions makes 
such a technique interesting for both applications and 
fundamental studies. The phenomenology observed for 
the junctions based on these GBs and Scanning SQUID 
Microscopy investigations demonstrate the absence of 7r- 
loops, as we expect from their microstructure. As a con- 
sequence higher values of the IcRn values, a Fraunhofer 
like dependence of Ic on the magnetic field and lower 
values of the low frequency flux noise, when compared 
with 45° c-axis tilt GBs, have been measured. These 
features are important tests to employ junctions for ap- 
plications. Scanning SQUID Microscopy investigations 
also gave evidence of "fractional" vortices in the pres- 
ence of impurities. Finally, we extended the biepit axial 
process to other types of GB by using different seed layers 
to obtain junction configurations where 7r loops can be 
controllably produced. We shall not dwell on conceptual 
principles and actual feasibility of qubit devices. Instead 
we discuss the importance of the biepitaxial technique in 
having "0" and w 7r" loops on the same chip. This makes 
the biepitaxial technique more versatile and promising 
for circuit design. 



of the GB interface, we have also successfully employed 
vicinal substrates. However, most of the transport prop- 
erties presented in this paper refer to samples not using 
vicinal substrates. Detailed structural investigations on 
these GBs, including Transmission Electron Microscopy 
(TEM) analyses, have been performed and the results 
have been presented elsewhere. 10,16 

Depending on the patterning of the seed layer and the 
YBCO thin film, different types of GBs ranging from the 
two ideal limiting cases of 45° a-axis tilt and 45° a-axis 
twist have been obtained (see Fig.l). The intermediate 
situation occurs when the junction interface is tilted at 
an angle a different from 0 or ir/2 with respect to the a- 
or b-axis of the [001] YBCO thin film. In all cases, the or- 
der parameter orientations do not produce an additional 
7r phase shift along our junction, in contrast with the 
45° asymmetric [001] tilt junctions. As a consequence, 
no 7T loops should occur independently of the details of 
the interface orientation. In Fig. 1 we consider ideal in- 
terfaces and neglect meandering of the GBs or interface 
anomalies that will be considered below. The CeC>2 seed 
layer may produce a more complicated GB structure, in 
which a 45° c-axis tilt accompanies the 45° a-axis tilt or 
twist (see Fig.2a). 15 In this case, as shown in Fig.2b, tt 
loops should occur in analogy with the traditional biepi- 
taxial junctions based on 45° c-axis tilt GBs. In both 
Figs. 1 and 2 we display the possible d x 2_ y a -wave or- 
der parameter symmetry in the junction configuration. 
Junctions were typically 4 microns wide. We also per- 
formed systematic measurements on SQUIDs based on 
the structure employing MgO as a seed layer and SrTiC>3 
as a substrate. DC SQUIDs in different configurations 
and with loop inductance typically ranging from 10 to 
100 pH have been investigated. The typical loop size 
leading to the 10(100) pH inductance is approximately 
lOVm 2 (lOVm 2 ). 



II. DEVICES: CONCEPTS AND FABRICATION 
PROCEDURE 

As mentioned above, the biepitaxial technique allows 
the fabrication of various GBs by growing different seed 
layers and using substrates with different orientations. 
We have used MgO, Ce02 and SrTi03 as seed layers. The 
MgO and Ce02 layers are deposited on (110) SrTi03 sub- 
strates, while SrTi03 layers are deposited on (110) MgO 
substrates; in all these cases the seed layers grow along 
the [110] direction. Ion milling is used to define the re- 
quired geometry of the seed layer and of the YBCO thin 
film respectively, by means of photoresist masks. YBCO 
films, typically 120nnv:n thickness, are deposited by in- 
verted cylindrical magnetron sputtering at a temperature 
of 780° C. YBCO grows along the [001] direction on MgO 
(substrates or seed layers) and on the Ce02 (seed layers), 
while it grows along the [103]/[013] direction on SrTi0 3 
(substrates or seed layers). In order to select the [103] or 
[T03] growth and to ensure a better structural uniformity 



III. EXPERIMENTAL RESULTS 

A. Biepitaxial junctions employing MgO seed layers 

In this section we attempt to cover most of the phe- 
nomenology of the transport properties of 45° a-axis tilt 
and twist biepitaxial junctions. In Fig. 3, current vs volt- 
age (I-V) characteristics of a typical biepitaxial junction 
are given for various temperatures close to the critical 
temperature. In the inset the corresponding I-V charac- 
teristic at T = 4.2 K is reported. They are closely de- 
scribed by the resistively-shunted-j unction (RSJ) model 
and no excess current is observed. Nominal critical cur- 
rent densities Jc of 5 x 10 2 A/cm 2 at T = 77 K, and of 9 x 
10 3 A/cm 2 at T = 4.2 K have been measured respectively. 
The Ryv value (3.2 Q) is roughly independent of the tem- 
perature for T < Tc, providing a normal state specific 
conductance a N - 70 (/xncm 2 )" 1 . The maximum work- 
ing temperature Tc of this device was 82 K. In this case 
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I c Rn is 1.3 mV at T = 4.2 K. These values typically 
ranged from 1 mV to 2 mV at T = 4.2 K. They are larger 
for the corresponding Jc values than those provided by 
conventional biepitaxials, and are of the same order of 
magnitude as in GB bicrystal and step edge junctions. 10 
While the values of critical current density and normal 
state specific conductance in the tilt case are quite dif- 
ferent from the twist case, the IcRn values are approxi- 
mately the same, for both. Moreover IcR-n does not scale 
with the critical current density. 10 In the tilt cases Jc « 
0.5-10 10 3 A/cm 2 and <tn « 1-10 (/iflcm 2 )" 1 are mea- 
sured at T = 4.2 K respectively. Twist GBs junctions 
are typically characterized by higher values of Jc in the 
range 0.1-4.0 x 10 5 A/cm 2 and of as in the range 20-120 
(/xHcm 2 )" 1 (at T = 4.2 K). For the twist case deviations 
from the RS J model are more marked as a result of higher 
critical current densities. For high values of Jc GB junc- 
tions do not present any clear modulation of the critical 
current as a function of the magnetic field. 

A demonstration of the possibility of tailoring the criti- 
cal current density and of the different transport regimes 
occurring in the tilt and twist cases has been given by 
measuring the properties of junctions with different ori- 
entations of the GB barrier on the same chip. By pattern- 
ing the seed layer as shown in Fig. 4a, we could measure 
the properties of a tilt junction and of junctions whose in- 
terface is tilted in plane by an angle a = 30°, 45° and 60° 
with respect to the a- or b-axis of the [001] YBCO thin 
film respectively. In all cases the order parameter orien- 
tations do not produce an additional 7r phase shift along 
our junction, in contrast with the 45° [001] tilt junctions, 
and no 7r loops should occur. We measured the expected 
increase of the critical current density with increasing an- 
gle, which corresponds on average to an increase of the 
twist current component. The values measured at T = 
4.2 K are reported in Fig. 4a and range from the min- 
imum value Jc = 3 x 10 2 A/cm 2 in the tilt case to the 
maximum Jc = 10 4 A/cm 2 corresponding to an angle of 
60°, for which the twist component is higher. The con- 
sistency of this result has been confirmed by the values of 
normal state resistances, which are higher in the tilt case 
and decrease with increasing a. The IcRtv values are 
about the same for all the junctions independently of the 
angle a. In Fig. 4b the I-V characteristics measured at 
T = 4.2 K, corresponding to the junctions of Fig.4a, are 
shown for approximately the same voltage range. Devia- 
tions from RSJ behavior appear for higher values of the 
critical current density (a = 60°). These results demon- 
strate that the grain boundary acts as a tunable bar- 
rier. This possibility of modifying the GB macroscopic 
interface plane by controlling the orientation of the seed 
layer's edge is somehow equivalent to the degree of free- 
dom offered by bicrystal technology to create symmetric 
or asymmetric GBs, with the advantage of placing all the 
junctions on the same substrate. The 45° a-axis tilt and 
twist GBs and the intermediate situations can represent 
ideal structures to investigate the junction physics in a 
wide range of configurations. The anisotropy of the (103) 



films and the possibility to select the orientation of the 
junction interface by suitably patterning the seed layer, 
and eventually the use of other seed layers which produce 
different YBCO in plane orientations, allow the fabrica- 
tion of different types of junctions and the investigation 
of different aspects of HTS junction phenomenology. In 
particular we refer to the possibility of changing the tun- 
neling matrix elements (by selecting the angle a) and 
to use the anisotropy of the layered structure of YBCO 
properties and of the order parameter symmetry. 

The study of the junction properties in the presence 
of an external magnetic field H is a fundamental tool 
for the investigation of the Josephson effect in the vari- 
ous junctions, as well as a test of junction quality. 17 We 
observe modulations of the critical current Ic following 
the usual Fraunhofer-like dependence. The Ic (H) pat- 
terns are mostly symmetric around zero magnetic field, 
and in all samples the absolute maximum of Ic occurs 
at H=0. The presence of the current maximum at zero 
magnetic field is consistent with the fact that in our junc- 
tion configuration the order parameter orientations do 
not produce an additional 7r phase shift, in contrast with 
the 45° [001] tilt GB junctions. 14 ' 10 Some examples are 
given in Fig. 5, where the magnetic pattern relative to 
a SQUID and a single junction at T = 4.2 K are shown 
respectively. In the former case we can also distinguish a 
smaller field modulation (with a period of 8 mG) which 
corresponds to the SQUID modulation (inset a). In the 
latter case the I-V characteristics are reported for differ- 
ent magnetic fields (inset b). Despite the Fraunhofer-like 
dependence, some deviations are evident, in agreement 
with most of the data available in literature. 

For sake of completeness we also acknowledge some 
work we carried out by investigating Fiske steps as a 
function of H in other junctions, giving some evidence of 
a dielectric-like behavior 18 of some of the layers at the 
junction interface. We already reported about this work 
elsewhere. 19 The Fiske steps do not depend on the use 
of a particular substrate, since they have been observed 
in junctions based both on SrTiOa and MgO substrates. 
Typical values of the ratio between the barrier thickness 
t and the relative dielectric constant e r range from 0.2 
nm to 0.7 nm. Considerations on the dependence of Ic 
on the temperature (T) can be also found in Ref. 19 . In 
junctions characterized by lower critical current densities, 
Ic tends to saturate at low temperatures, in contrast to 
those characterized by higher critical currents, for which 
there is a linear increase. 10,20 



B. Scanning SQUID microscopy on biepitaxial 
junctions with MgO seed layer 

Figure 6 is a scanning SQUID microscope 21 image of 
a 200x200/xm 2 area along a grain boundary separating a 
(100) region from a (103) region (as labelled in the fig- 
ure) of a thin YBCO biepitaxial film grown as described 
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above. The position of the grain boundary is indicated 
by the dashed line. The image was taken at 4.2K in liquid 
helium with an octagonal SQUID pickup loop 4 microns 
in diameter after cooling the sample in a few tenths of a 
jxT externally applied magnetic field normal to the plane 
of the sample. The grey-scaling in the image corresponds 
to a total variation of 0.13$o of flux through the SQUID 
pickup loop. Visible in this image are elongated inter- 
layer Josephson vortices in the (103) area to the right, 
and "fractional" vortices in the (100) area to the left, of 
the grain boundary. Fits to the interlayer vortices give 
a value for the c-axis penetration depth of about 4/jm. 
The "fractional" vortices are spontaneously generated in 
the (100) film, regardless of the value of external field 
applied. 22 Temperature dependent scanning SQUID mi- 
croscope imaging shows that this spontaneous magneti- 
zation, which appears to be associated with defects in the 
film, arises when the film becomes superconducting. 4,23 
Although it is difficult to assign precise values of total 
flux to the "fractional" vortices, since they are not well 
separated from each other, fits imply that they have less 
than $ 0 of total flux in them, an indication of broken 
time-reversal symmetry. Although there is apparently 
some flux generated in the grain boundary region, the 
fact that these SQUIDs have relatively low noise seems 
to indicate that this flux is well pinned at the tempera- 
tures at which the noise measurements were made. These' 
results are consistent with the absence of 7r loops along 
the grain boundary. 



C. Biepitaxial junctions employing Ce02 seed layers 



The CeC>2 seed layer, as anticipated in section II, may 
produce an artificial GB that can be seen as a result of 
two rotations: a 45° [100] tilt or twist followed by a 45° 
[001] tilt around the c-axis of the (001) film. For this 
junction configuration a d-wave order parameter sym- 
metry would produce 7r-loops, as shown in Fig. 2. We 
notice that such 7r-loops are structurally different from 
those usually obtained by the 45° [001] tilt GB junctions 
based on the traditional biepitaxial and bicrystal tech- 
niques. Due to the microstructure we expect especially 
in the [100] tilt case low critical current densities and 
high normal state resistances. We found that the depo- 
sition conditions to select the uniform growth of YBCO 
45° tilted around the c-axis of the (001) film are critical. 
Preliminary measurements realized on tilt- type junctions 
with a Ce02 seed layer gave evidence of Josephson cou- 
pling in these GBs. The measured IcR-vv values are from 
200 /xV to 750 jiV and are in the typical range of the GBs 
Josephson junctions. 



D. Biepitaxial SQUIDs employing MgO seed layers 



In this section we report on the characterization of dc- 
SQUIDs which are to our knowledge the first employing 
the GBs discussed above. 24 These SQUIDs exhibit very 
good properties, and noise levels which are among the 
lowest ever reported for biepitaxial junctions. 24 Apart 
from implications for applications, these performances 
are important for the study of the transport properties 
of HTS Josephson junctions. In Fig. 7 we show the mag- 
netic field dependence of the voltage at 77 K for different 
values of the bias current for a dc-SQUID with an induc- 
tance of 13 pH. At this temperature IqRn is about 20 
/zV. The corresponding value of the screening parameter 
/?=2LIc/$o is 0.03. In general low 0 values are manda- 
tory to avoid the influence of asymmetric inductances in 
SQUID properties, and this has been crucial for experi- 
ments designed to study the order parameter symmetry. 7 
The presented curves are quite typical. These SQUIDs 
usually work in a wide temperature range from low tem- 
peratures (4.2 K ) up to temperatures above 77 K. The 
maximum working temperature was in this case 82 K. 
The achieved magnetic flux-to- voltage transfer functions 
V* = dV/dQ, where V and $ are the voltage across the 
device and the applied magnetic flux in the SQUID loop 
respectively, are suitable for applications. For instance 
at T = 77 K an experimental value of the SQUID ampli- 
tude voltage modulation AV of 10.4 mV was measured, 
corresponding to V* = 36.9 /iV/$ 0 - 24 Steps of different 
nature have been recurrently observed in the I-V charac- 
teristics in the washer and hole configurations and char- 
acterized also in terms of the magnetic field dependence 
of the voltage at different values of the bias current. 

The noise spectral densities of the same dc-SQUID 
have been measured at T = 4.2 K and T = 77 K using 
standard flux-locked-loop modulated electronics. The en- 
ergy resolution e = S<j,/2L (with S* being the magnetic- 
flux-noise spectral density) at T = 4.2 K and T = 77 K is 
reported in Fig. 8. At T = 4.2 K and 10 kHz, a value of 
S* = 3 fi&o/VHz has been measured, corresponding to 
an energy resolution e = 1.6xl0~ 30 J/Hz. This value is 
the lowest reported in the literature for YBCO biepitaxial 
SQUIDs. Moreover, the low frequency 1/f flux noise spec- 
tral density at 1 Hz is more than one order of magnitude 
lower than the one reported for traditional biepitaxials, 
as is also evident from the comparison with data at T = 
4.2 K of Ref. 25 . The lower values of low frequency noise 
are consistent with the absence of 7r- loops on the scale of 
the faceting for these types of GBs, as clearly shown by 
Scanning SQUID Microscopy results. The 7r-loops pro- 
duce some types of spontaneous magnetic flux in the GB 
region, which among other effects tends to degrade the 
SQUID's noise levels. 14 
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IV. BIEPITAXIAL JUNCTIONS FOR 
EXPERIMENTS ON THE SYMMETRY OF THE 

ORDER PARAMETER AND FOR A 
DEVELOPMENT OF CONCEPTS FOR QUBITS 



The particular junction configurations investigated in 
this work allow some consideration of the possible impact 
of these types of junctions on the study of the Josephson 
effect and the order parameter symmetry in YBCO and 
on the development of concepts for devices. 1,6t9 7 We first 
recall that the biepitaxial technique can provide circuits 
composed completely of junctions without any 7r-loops 
(see Fig. 9a). By varying the interface orientation with 
respect to the [103] electrode orientation, the junction 
properties can be adjusted. On the other hand the tra- 
ditional biepitaxial technique, 11 producing 45° [001] tilt 
GBs (see Fig. 9b) or the types of junctions described in 
the previous section by using Ce02 (see Fig. 9c), can 
controllably generate 7r-loops on macroscopic scales. In 
these schemes we use a corner geometry with a 90° an- 
gle. This angle a can be obviously tuned to enhance 
the effects related to the phase shift (see dashed line in 
Fig.9b) and this change is particularly easy to realize by 
using the biepitaxial technique. 

In this section we focus our attention mainly on the 
feasibility of the biepitaxial junctions to obtain the dou- 
bly degenerate state required for a qubit. In Ref. la the 
design is based on quenching the lowest order coupling 
by arranging a junction with its normal aligned with 
the node of the d-wave order parameter, thus produc- 
ing a double periodic current-phase relation. It has been 
shown that the use of 7r phase shifts in a superconducting 
phase qubit provides a naturally bistable device and does 
not require external bias currents and magnetic fields. 16 
The direct consequence is the quietness of the device 
over other designs. A n junction provides the required 
doubly degenerate fundamental state, which also mani- 
fests itself in a doubly periodic function of the critical 
current density as a function of the phase. 8 The same 
principle has been used in small inductance five junction 
loop frustrated by a 7r-phase shift. 16 This design pro- 
vides a perfectly degenerate two-level system and offers 
some advantages in terms of fabrication ease and perfor- 
mance. HTS may represent a natural solution for the 
realization of the required 7r-phase shift due to the pair- 
ing symmetry of the order parameter and, therefore, due 
to the possibility of producing 7r phase shifts. Experi- 
mental evidence of YBCO 7r -SQUIDs has been given by 
employing the bicrystal technique on special tetracrystal 
substrates. 7 The biepitaxial technique, beyond provid- 
ing junctions with opportune properties, would guaran- 
tee the versatility necessary for the implementation of a 
real device, as shown below. As a matter of fact, we no- 
tice that our technique allows the realization of circuits 
where 7r- loops can be controllably located in part of the 
substrate and separated from the rest of the circuit based 
on "0" -loops, i.e. junctions where no additional n phase 



shifts arise. This can be easily made by depositing the 
MgO and "Ce02 seed layers on different parts of the sub- 
strate, which will be also partly not covered by any seed 
layer. 

As a test to show how the biepitaxial junctions could 
be considered for preliminary tests and device implemen- 
tation for quantum computing without the topological 
restriction imposed by the bicrystal technique, we refer 
to the structures proposed in Ref. 1 as exemplary circuits. 

The former is composed by a s-wave (S)- d-wave (D)- s- 
wave (S') double junction connected with a capacitor and 
an ordinary "0" Josephson junction based on s-wave su- 
perconductors (the S-D'-S junction generates the doubly 
degenerate state). The latter consists of a five junction 
loop with a 7T junction. Our technique would combine 
the possibility of placing the ordinary "0" junctions cor- 
responding to the MgO seed layer and to exploit the pos- 
sible doubly degenerate state of asymmetric 45° GB junc- 
tions corresponding to the Ce02 seed layer to replace the 
S-D-S' system or the 7r junction respectively. Our struc- 
ture would be obviously composed only of HTS. In Figs. 
10a and 10b we show how devices for instance such as 
those proposed in Ref. 1 could be obtained by employing 
the biepitaxial technique respectively. The application to 
the five junction loop is straightforward (Fig. 10b) and the 
advantages of this structure have been already discussed 
in Ref. 16 . The biepitaxial technique can offer possible 
alternatives for the realization of the structures above. 
In particular the double junctions of the original S-D- 
S' system can be also replaced by a D'-D-D" structure 
(Fig. 10c) by exploiting our technique, in contrast to the 
bicrystal technology which could not give this possibil- 
ity. Such a configuration could offer some advantages, if 
we consider that asymmetric 45° bicrystal GB Joseph- 
son junctions did not give systematic evidence of the 
doubly degenerate state. The doubly degenerate state 
seems to occur only in high quality low transparency GB 
junctions 8,26 and it is known that S-I-D junctions do not 
have double periodicity of the critical current as a func- 
tion of the phase. 26 A consequence of a possible nodeless 
order parameter 4,23 at the D-D' GB interface could be 
a closer similarity with a S-I-D junction with loss of the 
doubly degenerate state. If this is the case, we spec- 
ulate that the double junctions structure for symmetry 
reasons would produce a leading term in the Josephson 
coupling energy of the form Ed cos20 (double periodic) 
and that the possible dipolar component of the magnetic 
field would be almost completely compensated in this 
configuration. 16 This can be considered as an attempt to 
construct a "microscopic" 20-junction. We finally notice 
that the topological advantages offered by the biepitax- 
ial junctions would therefore be crucial in both the cases 
considered for the realization of the structure in Fig. 10, 
and important to reduce de- coherence effects. Bicrystal 
substrates would in fact impose on the circuit additional 
junctions required by the circuit design and, as a con- 
sequence, generate additional noise and de-coherence in 
the device. 
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V. CONCLUSIONS 

The performance of the presented junctions and 
SQUIDs demonstrates that significant improvements in 
the biepitaxial technique are possible, and the resulting 
devices have potential for applications. We have pre- 
sented a phenomenology that is consistent with the ex- 
pected absence of 7r-loops in 45° [100] tilt and twist grain 
boundaries junctions. The use of a Ce02 rather than a 
MgO seed layer can produce 7r-loops in the same junction 
configurations. The versatility of the biepitaxial tech- 
nique has been recently used to obtain different types 
of grain boundaries. The advantage of placing junctions 
in arbitrary locations on the substrate without imposing 
any restrictions on the geometry, and the ease of obtain- 
ing different device configurations by suitably patterning 
the seed layer, make the biepitaxial technique competi- 
tive for the testing of new concept devices, such as those 
based on 7r-loops. Some simple examples of situations 
in which 7r-loops can be suitably produced in specific 
locations of a more complicated circuit have also been 
discussed. 
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FIG. 1. A schematic representation of the artificial grain 
boundary structure. The boundary is obtained at the inter- 
face between the [001] oriented YBCO film grown on the [110] 
MgO seed layer and the [103] YBCO film grown on the bare 
[110] STO substrate. In contrast with the 45° [001] tilt bicrys- 
tal junctions, in this case the order parameter orientations do 
not produce an additional n phase shift. 
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FIG. 2. The CeC>2 seed layer produces an artificial GB 
that can be seen as a result of two rotations: a 45° [100] tilt 
or twist followed by a 45° tilt around the c-axis of the (001) 
film. For this junction configuration a d-wave order parameter 
symmetry would produce 7r -loops. 

FIG. 3. Current vs voltage (I-V) characteristics of the 
biepitaxial junction for temperature close to the critical tem- 
perature. In the inset the I-V curve at T = 4.2 K is shown. 

FIG. 4. a) Scheme of the seed layer patterning, which al- 
lows the measuremment on the same chip of the properties 
of a tilt junction and of junctions whose interface is tilted in 
plane of an angle a = 30° , 45° and 60° with respect the a- or 
b-axis of the (001) YBCO thin film respectively, b) The I-V 
characteristics (measured at T = 4.2 K) of the microbridges 
reported in Fig. 4a. 

FIG. 5. Magnetic-field dependence of the critical current 
of a [100] tilt biepitaxial dc-SQUID. The absolute maximum 
is observed for zero field. A double-period modulation is ob- 
served. The longer period modulation is the diffraction pattern 
due to the magnetic field sensed by a single junction, while 
the shorter period SQUID modulation is shown more clearly 
in the inset (a). In the inset (b) I-V curves are shown as a 
function of an externally applied magnetic field at T = 4.2 K. 
A typical Fraunhofer-like dependence is evident. 



FIG. 6. Scanning SQUID microscope image of a 
200x200/ira 2 area along a grain boundary separating a (100) 
region from a (103) region of a thin YBCO biepitaxial film 
grown. The position of the grain boundary is indicated by 
the dashed line. 



FIG. 7. Magnetic field dependence of the voltage of a [100] 
tilt biepitaxial dc-SQUID at 77 K for different values of the 
bias current. 

FIG. 8. Magnetic flux noise spectral densities of a [100] 
tilt biepitaxial SQUID at T=77 K and T=4.2 K. The SQUID, 
with an inductance L=13 pH, was modulated with a standard 
flux-locked-loop electronics. The right axis shows the energy 
resolution. Data at T = 4.2 K are compared with results on 
SQUIDs based on [001] tilt biepitaxial junctions from Ref. 25. 

FIG. 9. a) 3-dimensional view of a SQUID based on 45° 
[100] tilt and twist GBs; no 7r-loops should occur, b) Top 
view of tt-SQUID based on 45° [001] tilt GBs. c) 3-dimen- 
sional view of a 7T-SQUID based on GBs resulting from two 
rotations: a 45° [100] tilt or twist followed by a 45° [001] tilt 

FIG. 10. Scheme of the qubit structure proposed in Ref.l 
designed using the biepitaxial grain boundaries proposed in 
the paper. The double junctions of the original S-D-S' system 
can be also replaced by D'-D-D". 
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The current-phase relation (CPR) for asymmetric 45° Josephson junctions between two d-wave 
superconductors has been predicted to exhibit an anomalous periodicity. We have used the single- 
junction interferometer to investigate the CPR for this kind of junctions in YBa2Cu307-x thin films. 
Half-fluxon periodicity has been experimentally found, providing a novel source of evidence for the 
d-wave symmetry of the pairing state of the cuprates. 
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There is growing evidence in favor of the d x 2_ y 2-wave 
symmetry of the pairing state of the high-temperature 
superconductors. 1 An unconventional pairing state re- 
quires the existence of zeros of the order parameter in 
certain directions in momentum space. Thermodynamic 
and spectroscopic measurements do indeed suggest their 
existence, but by themselves they do not exclude conven- 
tional s- wave pairing with nodes. 1 Direct evidence for 
the d-wave pairing state is provided by phase-sensitive 
experiments, which are based on the Josephson effect. 2 
Quite generally, the current-phase relationship (CPR) of 
a Josephson junction, I {if), is an odd periodic function 
of if with a period 2ir. 3 Therefore I(tp) can be expanded 
in a Fourier series 



I(tp) = I\ sin (p + I2 sin 2<p + . 



(i) 



In the tunnel limit we can restrict ourselves to the first 
two terms in Eq. (1). Since the order parameter is 
bound to the crystal lattice, I(<p) of a weak link de- 
pends on the orientation of the d-wave electrodes with 
respect to their boundary. The existing phase-sensitive 
experiments exploit possible sign changes of I\ between 
different geometries. 2 In this Letter we present a new 
phase-sensitive experimental test of the symmetry of the 
pairing state of the cuprates. Namely, in certain geome- 
tries, the I\ term should vanish by symmetry. In such 
cases, the CPR should exhibit an anomalous periodicity. 

Let us analyze the angular dependence of I\, 2 in a junc- 
tion between two macroscopically tetragonal d-wave su- 
perconductors. As emphasized in Ref. 4, also heavily 
twinned orthorhombic materials such as YBa2Cu307_ z 
belong to this class, if the twin boundaries have odd sym- 
metry. We consider first an ideally flat interface between 
the superconducting electrodes. Let 0\ {$2) denote the 
angle between the normal to the grain boundary and the 
a axis in the electrode 1 (2), see Fig. 1. If we keep only 
the lowest-order angular harmonics, the symmetry of the 
problem dictates that 4 

(2) 



em dictates mat' 
Ii = I c cos 20i cos 20 2 -r 1 3 sin 20i sin 26 2 . 



The coefficients I c ,Is are functions of the barrier 
strength, temperature T, etc. The I2 term results from 
higher-order tunneling processes and we neglect its weak 
angular dependence. It is seen from Eq. (2) that the cri- 
terion for the observation of an anomalous CPR, I\ — 0, 
is realized for an asymmetric 45° junction, i.e. a junction 
with 6\ = 45° and 0 2 — 0. For an interface which is not 
ideally flat, 0* — 0i{x) are functions of the coordinate x 
along the junction. Ji = 0 remains valid even in this 
case, if the average values (0i(x)) = 45° and (0 2 (x)) = 0. 
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FIG. 1. Sketch of the interferometer (not in scale). 

The I2 term is present also in weak links based on 
conventional s-wave superconductors but for all known 
types of weak links Ih/h \ < 1- For instance, for a tunnel 
junction \h/h\ ^ L For a SNS junction, J oc sin<p/2 
at T — 0, 5 and the Fourier expansion Eq. (1) leads to 
h/h ~ -2/5. Therefore, a possible experimental ob- 
servation of \I 2 /h\ ^ 1 in an asymmetric 45° junction 
provides direct evidence of d-wave symmetry of the pair- 
ing state in the cuprates. 

We have investigated the CPR of YBa 2 Cu 3 0 7 - I 
thin film bicrystals with asymmetric 45° [001]-tilt 
grain boundaries 6 " 8 as sketched in Fig. 1, using a 
single-junction interferometer configuration in which the 
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Josephson junction is inserted into a superconducting 
loop with a small inductance L. In a stationary state 
without fluctuations, the phase difference <p across the 
junction is controlled by applying external magnetic flux 
$ e penetrating the loop: <p - <p e - /?/(</?)• Here <p e = 
2n$ e /$o is the external flux normalized to the flux quan- 
tum $o = 2.07 x 10~ 15 Tm 2 . The CPR is written as 
j(<p) = Iof(<p), where Iq is the maximal Josephson cur- 
rent. 0 = 27tLIq/$q is the normalized critical current. 
In order to obtain the CPR for the complete phase range 
-7r < < 7r the condition 0 < 1 has to be fulfilled, 
because for 0 > 1 the curve tp{ip e ) becomes multivalued 
and there are jumps of <p and a hysteresis for a sweep 
of if e . Following Ref. 3, we express the effective induc- 
tance of the interferometer using the derivative /' with 
respect to <p as L int = L[l+l//'((p)]. The inductance can 
be probed by coupling the interferometer to a tank cir- 
cuit with inductance Lr, quality factor Q, and resonance 
frequency wo- 9 External flux in the interferometer is pro- 
duced by a current Idc + hf in the tank coil and can be 
expressed as y e = 27r(Jd c + J r /)M/$o = <A*c + ¥V/, where 
M 2 = k 2 LL,T, and A; is a coupling coefficient. Taking into 
account the quasiparticle current in the presence of volt- 
age V across the junction the phase difference is given 
by the relation <p = tpdc + <Prf - /?/(<?) - 27rr(cp)Vy$ 0 , 
where r(ip) = L/Rj(<p) and Rj(tp) is the resistance of 
the junction . In the small-signal limit <p r f 1 and in 
the adiabatic case ljt 1, keeping only the first-order 
terms, the effective inductance L e f / of the tank curcuit- 
interferometer system reads 

Thus the phase angle a between the driving current and 
the tank voltage U at the resonant frequency of the tank 
circuit u>o is 



tana(v?) = 



k 2 QPf'(<p) 
1 +/?/'(¥>)' 



(3) 



Using the relation [1 + /3/'(<^)]<fy> = d<pd c valid for (p r f < 
1 and wr < 1, one can find the CPR from Eq. (3) by 
numerical integration. 

The advantage of the measurement of the CPR of an 
asymmetric 45° junction with respect to the by-now stan- 
dard phase-sensitive tests of pairing symmetry based on 
the angular dependence of h is twofold. First, it avoids 
the complications of the analysis of experiments caused 
by the presence of the term I B * Second, a flux trapped 
in the SQUID does not invalidate the conclusions about 
the ratio \h/h\ and hence about the pairing symmetry, 
while this is not the caso in standard phase-sensitive tests 
of the d-wave symmetry of the pairing state. 10 

The films of thickness 100 nm were fabricated us- 
ing standard pulsed laser deposition on (001) oriented 
SrTiCb bicrystalline substrates with asymmetric [001] tilt 
misorientation angles 45° ±1°. They were subsequently 
patterned by Ar ion-beam etching into 4x4 mm 2 square 



washer single-j unction interferometer structures (Fig. 1). 
The widths of the junctions were 1 -5- 2 jim. The washer 
square holes had a side-length of 50 /im. This geometry 
of the interferometer gives L % 80 pH. The resistance of 
the junction is higher than 1 and the condition for the 
adiabatic limit ljt < 1 is satisfied. For measurements of 
a(<£d C ), several tank circuits with inductances 0.2 ~ 0.8 
j/H and resonance frequencies 16 -5- 35 MHz have been 
used. The unloaded quality factor of the tank circuits 
70 < Q < 150 has been measured at various tempera- 
tures. The coupling factor k was determined from the 
period AI dc of a{I dc ) using MAI dc = $o- Its value var- 
ied between 0.03 and 0.09. The amplitude of I r f was 
set to produce the flux in the interferometer lower than 
0.1$ 0 . 

The measurements have been performed in a gas-flow 
cryostat with a five-layer magnetic shielding in the tem- 
perature range 4.2 < T < 90 K. The experimental setup 
was calibrated by measuring interferometers of the same 
size with 24° and 36° grain boundaries. We have stud- 
ied 5 samples, out of which sample No. 1 exhibited the 
most anomalous behavior. Samples Nos. 2,3 were less 
anomalous and the remaining two samples had high crit- 
ical currents and their I(<p) was conventional. In Fig. 2 
we plot the phase angle a as a function of the dc current 
Ida for samples Nos. 1,2. The data for the 36° junction 
is shown for comparison. Note that at T = 40 K the 
periodicity of a(Idc) is the same for all samples. 
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FIG. 2. Left panel: Phase angle between the driving current 
and the output voltage measured for the sample No. 1 at dif- 
ferent temperatures as a function of the dc current Idc- The 
curves are shifted along the y axis and the data for T = 30 
and 40 K are multiplied by factor 4 for clarity. From top to 
bottom, the data correspond to T =4.2, 10, 15, 20, 30 and 40 
K. The data measured on 36° bicrystals ($i a; 36°, 6 2 ^ 0) 
at T =40 K in the same washer geometry are shown for com- 
parison (open circles). Right panel: The same for the sample 
No. 2. From top to bottom, the data correspond to T =4.2, 
10, 15 and 20 K. 

We assume that the period of a{I dc ) at T =40 K, 
A/dc =9-6 /iA, corresponds to A(f d c = 27r. In order 
to determine the CPR we take <pdc = 0 at a maximum 
or minimum of a. This is necessary in order to satisfy 
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= 0) = 0, as required by general principles. 3 The ex- 
perimentally observed shift of the first local extreme of 
a(Idc) from Idc = 0 (Fig. 2) can be due to flux trapped 
in the interferometer washer. Most probably, this flux re- 
sides in the long junction of the interferometer. The long 
junction does not play an active role because the Joseph- 
son penetration depth is much shorter than its length, 
and external fields producing Idc are smaller than its crit- 
ical field. Nevertheless the long junction sets the phase 
difference for Idc = 0 at the small junction. 
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Our experimental results can be understood as follows. 
Deviations from ideal geometry of the asymmetric 45° 
junction, {8\) = 45° + a.\ and (#2) = c*2, lead to a fi- 
nite value of I\. Thus, imperfections of the junction 
increase its critical current. For this reason we believe 
that samples Nos. 2-5 contain imperfections and from 
now we concentrate on nearly ideal junctions (such as 
sample No. 1) with |qi|,|q2| 1- For such junctions, 
the ratio h/h exhibits the following temperature de- 
pendence. For T — ► 0, \hlh \ > 1- The region T ~ T c 
can be analyzed quite generally within Ginzburg-Landau 
theory. Let the electrodes be described by (macro- 
scopic) order parameters Ai (2 = lAle'^ 12 . Then the 
phase-dependent part of the energy of the junction is 
E = a[AxA$ + H.C.] + 6[(A X A£) 2 + H.C.] + . . . where 
a, 6, . . . depend weakly on T. 11 Thus for T close to T c we 
estimate h oc |A| 2 oc (T c -T) and I 2 oc |A| 4 a (T c -T) 2 , 
leading to I/2/AI ^ 1- These expectations are qualita- 
tively consistent with the experimental data shown in 
Fig. 4. 
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FIG. 3. a) Josephson current through the junction for the 
sample No. 1 as a function of the phase difference <p, deter- 
mined from the data in Fig. 2. The statistics of a (<p) was 
improved by folding the data back to the interval (0, 7r) and 
taking an average. The symmetry cu(<^) = a(— (p) was as- 
sumed, b) The same for the sample No. 2. The data for the 
asymmetric 36° bicrystal at T = 40 K (open circles) is also 
shown. 

In Fig. 3, we show the CPR determined from the data 
in Fig. 2. For all curves we have performed a minimal 
necessary shift consistent with I(<p — 0) = 0. Thus, for 
the samples Nos. 1,2 we have assumed that at tpdc = 0 
a minimum of a(<pdc) is realized. For an interferometer 
with a conventional s-wave weak link (and also for the 36° 
junction), at tpdc — 0 a maximum of a(<pd c ) is realized. 
Note that the maximum (minimum) of a(ipdc) at <pdc = 0 
implies a diamagnetic (paramagnetic) response of the in- 
terferometer in the limit of small applied fields. In Fig. 4 
we plot the coefficients I\ and h determined by Fourier 
analysis of the CPR for the sample No. 1 at various tem- 
peratures. With decreasing T, I/2I grows monotonically 
down to T = 4.2 K, while the I\ component exhibits only 
a weak temperature dependence. 



FIG. 4. Temperature dependence of the Fourier expansion 
coefficients 7i,2 determined from the experimental data in 
Fig. 3a. Solid lines are the Fourier expansion coefficients for 
the numerical data in Fig. 5. 

So far, our discussion was based solely on symmetry 
arguments. Let us attempt a more quantitative analy- 
sis of our data now. Two different microscopic pictures 
of asymmetric 45° Josephson junctions between d-wave 
superconductors have been considered in the literature. 
The first picture assumes a microscopically tetragonal 
material and an ideally flat interface. 11 " 13 Within this 
picture, there are two contributions to the Josephson 
current. The first is due to bulk states and in the tun- 
nel limit it is well described by the Sigrist-Rice term 
I c in Eq. (2). 14 The second is due to mid-gap states 
which develop close to the surfaces of unconventional 
superconductors. 15 J(y?) for the sample No. 1 calculated 
according to the model of Ref. 12 is shown in Fig. 5. 
The experimental data can be fitted by a relatively broad 
range of barrier heights. However, if we require the 36° 
junction to be fitted by the same (or smaller) barrier 
height as for the 45° junction, we conclude the barrier 
must be rather low. 16 The T dependence of I(tp) requires 
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a choice of T c « 60 K in the non-selfconsistent theory 
of Ref. 12. The reduction from the bulk T c — 90 K is 
probably due to a combined effect of surface degradation 
and order-parameter suppression at the sample surface. 
The temperature dependence of the ratio of the 7r and 
27r periodic components in I(ip) is seen to be in qualita- 
tive agreement with experimental data in Fig. 3a. This 
is explicitly demonstrated in Fig. 4 where we compare 
the experimentally obtained ha with the results of the 
Fourier analysis of the curves in Fig. 5. The divergence of 
J 2 as T 0 is an artifact of the ideal junction geometry 
assumed in Ref. 12. If the finite roughness of the inter- 
face is taken into account, this divergence is cut off and 
the experimental data in Fig. 4 do indeed resemble the- 
oretical predictions for a rough interface. 13 However, the 
nonselfconsistent theory of Ref. 12 is unable to explain 
the experimentally observed steep CPR close to the min- 
ima of the junction energy (see Fig. 3a). In the limit of 
vanishing barrier height, the theoretical CPR does have 
steep portions, but these are located close to the maxima 
of the junction energy (see also Ref. 11). 



=3 



' 1 ' 

4.2 It /"\ 

10 K. / \ 

- - - -15 K / \ 

. 20 K / \ 

30 K / ^ - v \ 

*0 K . \ V 









-1.0 -0.5 



0.0 0.5 

(p/n 



1.0 



FIG. 5. I((p) calculated according to Eq. (64) of Ref. 12 for 
a junction with &i = 45.5°, 6 2 = 0, Ad = 1.5, k — 0.5, and 
T c = 60 K. I(<p) at T = 40 K for the 36° bicrystal (open 
circles) was calculated for the same parameters except for 
0i = 36°. 

In a different approach to the asymmetric 45° junction, 
one assumes a heavily twinned orthorhombic material 
(which is macroscopically tetragonal, however) and/or 
a meandering interface with $i = 9i(x). 17AS Hence the 
critical current density j c (x) is a random function with 
a typical amplitude (|j c (s)|) ~ Jc If the average criti- 
cal current along the junction (j c ) < j c > a spontaneous 
flux is generated in the junction, and I/2/J1I > 1. 17,18 In 
particular, for ($1) = 45° and (6 2 ) = 0, there- is an equal 
amount of parts having positive and negative j c , leading 
to (j c ) = 0 and I\ = OrNote that also within the picture 
of Refs. 17,18, the d-wave symmetry of the pairing state 
is crucial, otherwise the condition (j c ) ^ jc is difficult to 
satisfy. 



Our present understanding of 7(<p) in the asymmetric 
45° junction is only qualitative. There is considerable 
experimental evidence 6-8 that the grain boundary junc- 
tions are at most piecewise flat. However, we cannot say 
whether the shape of I(ip) is dominated by the mid-gap 
states in the microscopically flat regions, or by sponta- 
neous flux generation due to the spatial inhomogeneity 
of the junction. This issue requires further study. 

In conclusion, we have measured the magnetic field re- 
sponse of a single-junction interferometer based on asym- 
metric 45° grain-boundary junctions in YBa2Cu307- x 
thin films. Half-fluxon periodicity has been experimen- 
tally found, in agreement with theoretical predictions for 
d x 2_ y a-wave superconductors. Hence, our results provide 
a novel source of evidence for the d-wave symmetry of the 
pairing state in the cuprates. 
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Bi-epitaxial grain boundary junctions in YBa 2 Cu 3 0 7 
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We have developed a new way of making grain boundary junctions in YBa 2 Cu 3 0 7 thin films 
by controlling the in-plane epitaxy of the deposited film using seed and buffer layers. 
We produce 45° grain boundaries along photolithographically defined lines. The typical value 
of the critical current density of the junctions is 10M0 4 A/cm 2 at 4.2 K and 10 2 -10 3 A/ 
cm 2 at 77 K, while the rest of the film has a critical current density of 1-3 X 10 6 A/cm 2 at 77 
K. The current-voltage characteristics of the junctions show resistively shunted junction 
behavior and we have used them to fabricate dc superconducting quantum interference devices 
(SQUIDs) which show modulation at temperatures well above 77 K. This is the first 
planar high T c Josephson junction technology that appears readily extendable to high T e 
integrated circuits. 



Most microelectronic applications of the high T c su- 
perconductors will rely upon microbridges such as super- 
conductor-normal metal-superconductor (S-N-S) junctions 
or other weak link geometries instead of superconductor- 
insulator-superconductor (S-I-S) tunnel junctions. A num- 
ber of weak-link structures provide reliable critical current 
reduction in YBa 2 Cu 3 0 7 ( YBCO) thin films; grain bound- 
ary junctions. 1 " 3 traditional edge junctions, 4 microbridges' 
grown across a sharp substrate step, 5 and S-N-S-type junc- 
tions with Au, 6 PrBa 2 Cu 3 0 7 , 7 - 8 or SrTi, _ x Nb x 0 3 (Ref. 9) 
normal layers. Of these various weak-link geometries, grain 
boundary junctions work well at temperatures close to 
T a have reasonably high IcR n products, and show behavior 
explained by a resistively shunted junction model with a 
relatively uniform current density. To date, the high angle 
grain boundary junctions that have been used were formed 
either by fusing differently oriented SrTi0 3 substrates 1 or 
by patterning the randomly occurring grain boundaries in 
granular films. 2 * 3 Superconducting quantum interference 
devices (SQUIDs) made from the grain boundaries on Sr- 
TiOj bicrystals had reasonable signal and low noise 10 up to 
temperatures very close to the T c o{ the film. However, this 
SrTi0 3 bicrystal technique has a major drawback: it cannot 
be readily extended to integrated circuits. On the other 
hand, SQUIDs made from granular films are plagued by 
low yield and the presence of grain boundaries in the 
SQUID loop itself. These lead to excessive flux noise and 
hysteresis in the voltage-flux response of the SQUID. 

In this letter, we report for the first time the successful 
fabrication of high T c weak-link Josephson junctions oper- 
ating at and above 77 K that are made from grain bound- 
ary junctions whose locations are determined by conven- 
tional photolithography. This result demonstrates the 
feasibility of an extendable technology for fabricating high 
T c Josephson devices which avoids the problems associated 
with the devices made from granular thin films and on 
bicrystal substrates. We reproducibly fabricate 45° grain 
boundaries in YBCO films by controlling their in-plane 
epitaxy using seed and buffer layers deposited on /--plane 
sapphire substrates, a method we call "bi-epitaxy." The 
method is quite general, and can be used on a wide variety 
of substrates. 



Recently, we reported the existence of two in-plane 
epitaxial orientations of c-axis oriented epitaxial YBCO 
thin films on yttria-stabilized zirconia (YSZ), resulting 
from the poor lattice match between YBCO and YSZ. This 
observation led us to realize that it would be possible to 
make bi-epitaxial grain boundaries on such a substrate us- 
ing seed and buffer layers, or, in an alternative geometry, 
on a well lattice-matched substrate using poorly lattice- 
matched seed layers. For example, we previously re- 
ported 12 the growth of YBCO thin films on r-plane sap- 
phire using buffer layers such as SrTi0 3 , CaTi0 3 , YSZ, and 
MgO. We emphasized the importance of the in-plane epi- 
taxy of YBCO, particularly because high-angle grain 
boundaries in YBCO films behave as weak links, leading to 
low critical current densities and high surface resis- 
tances. 12,13 Our efforts towards minimizing the occurrence 
of th^se high-angle grain boundaries taught us how to con- 
trol them and led us to design and demonstrate a variety of 
45° weak-link junctions. One example takes advantage of 
the differing epitaxial directions of two different layers, 
MgO and SrTi0 3 , when they are deposited on an r-plane 
sapphire substrate. Using these two layers as a seed and a 
buffer layer, we have succeeded in fabricating bi-epitaxial 
45° grain boundaries in YBCO on sapphire substrates by 
the following process. 

Schematic views of our device structure are shown in 
Fig. 1. We first deposit 3-30 nm of epitaxial MgO as a seed 
layer on an r-plane sapphire substrate. MgO is produced by 
laser ablating a Mg metallic target in a 2-20 mTorr oxygen 
atmosphere with a substrate temperature of 300-700 °C. 
We then mask the MgO with conventional photoresist and 
remove it from, for example, half of the substrate by either 
Ar ion beam milling or chemical wet etching with dilute 
phosphoric acid. We then grow 10-100 nm of epitaxial 
SrTi6 3 buffer layer by laser ablation on both the exposed 
sapphire surface and on the patterned MgO seed layer. The 
deposition conditions are an oxygen pressure of 100-200 
mTorr and a substrate temperature of 7 10-760 °C. The 
SrTi0 3 film grows in two different orientations separated 
by a 45° grain boundary. In the growth plane the epitaxial 
relations are SrTi0 3 [ 1 10]//Al 2 O 3 [l 1 20] and 
SrTiO 3 [100]//MgO[100]//A] 2 O 3 [1120], as illustrated in 
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FIG. 1. Schematic view of the device structure with an MgO seed layer 
and a SrTiOj buffer layer on an r-piane sapphire substrate, (a) top view, 
(b) side view. 



Fig. 1. We then immediately deposit YBCO, which grows 
epitaxially everywhere and thereby reproduces the grain 
boundary in the SrTiOj buffer layer. Finally the YBCO is 
patterned into an appropriate geometry. 

The virtue of this structure is that, apart from the 
predefined 45° grain boundaries, there are no other high- 
angle grain boundaries in the YBCO film. It has been re- 
ported 12 that YBCO films grown on SrTiOj buffer layers 
on /--plane sapphire substrates have high critical current 
density and low surface resistance. In addition, the SrTiOj 
layer grown on the MgO buffer layer does not have any 
high-angle grain boundaries, as demonstrated in a ^ scan of 
the SrTiOj (101), although it is found 13 that YBCO de- 
posited on MgO may have a number of high-angle grain 
boundaries. The probable explanation for this effect is that 
the energy needed to nucleate 45* misoriented SrTiOj 
grains on MgO is much higher than in the case of YBCO 
grown on MgO. Hence the YBCO film on SrTiOj on MgO 
also lacks any high-angle grain boundaries and has a high 
critical current density. Microbridges formed either side of 
the predetermined grain boundary have critical current 
densities of 1-3 X 10 6 A/cm J at 77 K. In order to show that 
the critical current reduction is indeed due to the grain 
boundary rather than simply the step at the edge, we 
etched only part of the MgO, reducing its thickness on one 
side from 20 to 10 nm, such that YBCO/SrTi0 3 /MgO 
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FIG. 2. X-ray 6 scan of YBCO ( 103) of the unpatterned structure. The 
peaks labeled with "A w belong to YBCO/SrTKVAUO, and the peaks 
labeled with "B" come from YBCO/$rTiO»/MgO/AUO». 



remained on the sapphire substrate on both sides of the 
step.The critical density across this step remained larger 
than 10 6 A/cm 2 at 77 K. 

A 4> scan of the YBCO ( 103) before the last patterning 
process is shown in Fig. 2. The peaks occurring 45° apart 
indicate that half of the YBCO is rotated 45° relative to the 
other half. Suitable etching experiments have shown that 
the peaks labeled with "A" belong to YBCO/SrTiO/ 
A1 2 0 3 and the peaks labeled with "B" come from YBCO/ 
SrTiOj/MgO/Al 2 Oj. Cross-sectional transmission electron 
microscope images and microdiffraction patterns show that 
the 45* grain boundary begins at the end of the MgO seed 
layer. 14 

As the first example of a device using multiple junc- 
tions made with our process we fabricated a square washer 
SQUID 15 as shown in the inset of Fig. 3. It was patterned 
across the grain boundary by photolithography and wet 
etched with dilute phosphoric acid. The current-voltage 
characteristic (I-VC) of the SQUID at 4.2 K, which ex- 
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FIG. 3. Current (vertical )-voltage (horizontal) characteristics of a dc 
SQUID at 4.2 K. The inset shows the geometry of the SQUID. 
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FIG. 4. Modulation voltage of the current-biased dc SQUID vs applied 
magnetic flux at various temperatures, 

hibits resistively shunted Josephson junction behavior, is 
shown in Fig. . 3. Each junction was 4 fim wide and 0.2 fim 
thick. The critical current of about 10 /xA at 4.2 K trans- 
lates to about 10 3 A/cm 2 critical current density across the 
junctions. The /<J?„ product of the junctions is about 420 
pV at 4.2 K. The I-VC of the SQUID at 77 K shows no 
zero resistance part, although it is nonlinear. The small 
junction coupling energy I&t/lir compared to the thermal 
energy k B T may lead to thermally activated phase slippage 
across the junctions, resulting in resistance at all cur- 
rents. 16 Detailed transport properties of the junctions in- 
cluding the magnetic field dependence of the critical cur- 
rent will be reported elsewhere. 17 

The voltage across the current-biased dc SQUID as a 
function of applied field is shown for various temperatures 
in Fig. 4. The data were taken in a bandwidth of 0-10 Hz 
without any signal averaging. The observed period corre- 
sponds to the expected value from the geometry of the 
SQUID. We believe that the reduction of the modulation 
voltage at high temperatures jis due to the decreasing crit- 
ical current and the large inductance of the SQUID. We 
have observed dc SQUID operation at temperatures as 
high as 88 K and have operated such a SQUID in j flux- 
locked loop up to 83 K. The noise of this SQUID is similar 
to that reported in dc SQUIDs made using bicrystal sub- 
strates. 10 Detailed performance of the SQUID will be pub- 
lished elsewhere. 18 

A further example of small-scale integration using 
biepitaxial grain boundary junctions is provided by a flux 



shuttle 19 we have fabricated. It involves two dc SQUIDs as 
sensors and eleven more bi-epitaxial junctions as switching 
elements. Our work on such integrated circuits is continu- 
ing and will be reported later. 

In summary, we have reported the development of bi- 
epitaxial grain boundary junctions in YBa 2 Cu 3 0 7 at mul- 
tiple and predetermined locations by using only standard 
photolithographic techniques to control in-plane epitaxy 
with seed and buffer layers. The junctions have current- 
voltage characteristics that are well described by the resis- 
tively shunted junction model and dc SQUIDs fabricated 
using these junctions show modulation at temperatures as 
high as 88 K. By increasing the critical current density of 
the junctions and lowering the inductance of the SQUID, 
the performance of the SQUIDs should improve further, 
especially at 77 K. We believe that this technology, based 
on bi-epitaxial grain boundaries, can be extended to inte- 
grated circuits in the near further, and we are planning to 
demonstrate such circuits. 
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